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The phenomenon of diquark condensation at non-zero baryon density and zero temperature is analyzed in the 
strong coupling limit of lattice QCD. The results indicate that there is attraction in the quark-quark channel also 
at strong coupling, and that the attraction is more effective at high baryon density, but for infinite coupling it 
is not enough to produce diquark condensation. It is argued that the absence of diquark condensation is not a 
peculiarity of the strong coupling limit, but persists at sufficiently large finite couplings. 



A weak coupling analysis of QCD, both at the 
level of one gluon exchange as well as using the 
instanton induced interaction, reveals that the 
quark-quark scattering amplitude is attractive in 
the color anti-triplet channel and repulsive in the 
symmetric color sextet channel. The existence of 
such attractive force between quarks in the color 
anti-triplet channel is also supported by experi- 
ments, since hadron phenomenology is best ex- 
plained by models that assume baryons consiting 
of a quark-quark anti-triplet color bound state 
whose color is neutralized by the remaining quark 
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In the deconfined phase, at high baryon den- 
sities and low temperatures, the effect may be 
more dramatic. If the weakly interacting quarks 
tend to form a Fermi surface, the quark-quark 
attraction will render it unstable via the Cooper 
mechanism. Hence, the phenomenon of color su- 
perconductivity will take place: a diquark con- 
densate will appear in the ground state, a gap 
will be opened in the spectrum and the Fermi 
surface will be destroyed. Many exotic phenom- 
ena at high baryon density may happen as a con- 
sequence of diquark condensation, as color flavor 
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locking and unlocking, crystalline color supercon- 
ductivity, etc. P]. 

The diquark operator is not gauge invariant, 
and therefore general arguments forbide diquark 
condensation in the naive way: it must be un- 
derstood not as the spontaneous breaking of the 
gauge symmetry, but as a kind of Higgs mech- 
anism whose physical consequences are the phe- 
nomena we name color superconductivity. 

The arguments leading to the existence of di- 
quark condensation at high baryon density are 
based on weak coupling analysis -although some 
non perturbative effects are taken into account 
via instantons-, and thus are only valid at as- 
symptotically large densities. In addition, the 
study is performed using effective four quark the- 
ories which are not gauge invariance, and then it 
is difficult to trace the fate of gauge invariance. 
It is of utmost importance to study diquark con- 
densation beyond the weak coupling region, be- 
cause of the lack of reliable tools, since Monte 
Carlo simulations are not feasible at high density 
and low temperature due to the sign problem - 
however, see Hence, it is interesting to turn 
to the opposite regime, the strong coupling do- 
main. Recent attempts to study finite density 
QCD in the strong coupling limit using hamilto- 
nian techniques have been developed in 0] . Here, 
we present the results concerning diquark con- 
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densation obtained in [S] by using path integral 
techniques. 

Let us consider QCD -the gauge group is 
SU(3)- at zero temperature and finite baryon 
chemical potential, regularized on a euclidean 
four dimensional lattice with four flavor of stag- 
gered fermions. Let us define the meson and 
baryon fields 

M{x) ^ ^P''{x)i;''{x) , (1) 

B{x) = ^e,,,r{x)ij''{x)^P%x), (2) 
6 

B{x) = l:eabcr{x)^p''{x)r{x), (3) 
6 

and the diquark field 

D^ix) = eabc^\x)i^'{x), (4) 
D'^ix) = eabc^\x)r{x). (5) 

where summation over repeated color indexes 
a, b, c is understood. Note that B{x) = 
l/6'ip''{x)D''{x) and B{x) = l/6iJj''{x)D''{x). At 
strong coupling, the gauge field can be integrated 
out and what remains is a local action for the 
fermions, which, by gauge invariance, depends 
only on M{x),B{x), and B{x). The action is 
hence multilinear in the quark variables ip°'{x) 
and ■4>'^{x). 

The action can be converted in a bilinear in the 
quark fields at the expense of introducing a real 
colorless bosonic field, (p{x) and a color triplet 
complex bosonic field, 4>°'{x), through a series of 
Hubbard- Stratonovich transformations After 
integrating out the quark field, the lattice grand 
canonical partition function is given by an effec- 
tive bosonic action which depends only on the 
fields (p{x), (p"'{x) and (l)\{x), and which is invari- 
ant under SU(3) color gauge transformations, so 
that it actually depends only on the gauge invari- 
ant combination cl)\{x)(j)°'{x), and on 1^9(2;). 

The field 1^9(2;) indeed represents the meson field 
M{x) since, by introducing suitable sources in the 
action, it is easy to see that (p{x) = cM{x), where 
c is a constant. Thus, 1^9(2;) is an order parame- 
ter of chiral symmetry and its expectation value 
is proportional to the chiral condensate. Analo- 
gously, it can be seen that the field (p^ (x) -which 
is a color anti-triplet- is linked to the diquark field 



by D°'{x) = F{x) (t>a{x), where F{x) is a gauge in- 
variant functional of (p{x), (j)'^{x) and (j>]^{x). As 
expected, diquark condensation takes place if the 
expectation value of 4i°'{x) is non-zero. 

Thus, the field is a kind of Higgs field 

which is useful to study color superconductivity in 
the same way as the theory of Ginzburg-Landau 
describes electric superconductivity. Here, how- 
ever, the effective action for has been ob- 
tained exactly from QCD, albeit in the strong 
coupling regime, and there is no gauge field cou- 
pled to it. 

To study the phase diagram of the theory we 
used the semi-classical approximation. The ef- 
fective potential is given by the bosonic effec- 
tive action evaluated at constant fields 1^9(2;) = (p 
and |'/>(2;)p = v"^. Its absolute minimum de- 
scribes the equilibrium state at a given tempera- 
ture and chemical potential. Other local minima, 
if present, are associated with metastable states. 
In the chiral limit and at zero temperature we 
get the following minima of the classical effective 
potential 5 : 

v = 0, = ± (V33- 5)^^^ Va* 

w=V^, (^ = Vn (6) 

1,^0, ip^O V/i> 0.4416 

Note that the positions of the minima are inde- 
pendent of /i. 

There is a first order phase transition at fif, « 
1.557. For fi < the -degenerate owing to chi- 
ral symmetry- absolute minima are at d = and 
(f = ±(V33 — 5)^/^, corresponding to a phase 
with chiral symmetry spontaneously broken. For 
fi > the absolute minimum is at u = and 
(p = 0; chiral symmetry is restored in this phase. 
Obviously, the transition is first order. For ^ < fic 
the baryon density is zero, and three for fi > iic- 
Three quarks per lattice site is the maximum al- 
lowed by Pauli's principle, so that the transition 
separates a phase of zero baryon density from a 
phase saturated of quarks. The system at any in- 
termediate density is thermodynamically unsta- 
ble and splits into domains of zero and saturated 
baryon density. This behavior is an artifact of 
the strong coupling limit and has been observed 
in other strong coupling analysis at zero temper- 
ature and finite chemical potential. 
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Figure 1. Baryon density as a function of the 
chemical potential in the metastable state with 
diquark condensation. 

The minimum at v = \/5 and — is 
metastable in both phases. It describes a chiral 
symmetric state with a diquark condensate. The 
baryon density that correspond to this metastable 
state is a smooth function of the chiral conden- 
sate (see figure^ that interpolates between zero 
density and saturation. This means that at any 
baryon density a state with diquark condensation 
can be formed and have some short life until the 
system splits into its zero density and saturated 
domains. The presence of the metastable state 
signals the attraction in the quark-quark channel. 
At strong coupling such attraction is not strong 
enough to form a stable diquark condensate. The 
energy difference between this metastable state 
and the equilibrium state decreases with /i, indi- 
cating that the quark-quark attraction becomes 
the more effective the higher the baryon density 
(see figure 121. However, the energy difference is 
nonzero even at high fi (in the saturated phase). 

The crucial question is whether this metastable 
state becomes stable at sufficiently large density 
at finite coupling. A partial answer is the follow- 
ing: assuming that the strong coupling expan- 
sion gives meaningful results, diquark condensa- 
tion cannot take place at sufficiently large cou- 
plings, g, whatever the chemical potential. The 
reason is that, as we have seen, the energy differ- 
ence between the state with diquark condensation 
and the stable state remains positive for any value 
of the chemical potential in the strong coupling 
limit. A correction to the effective potential of 



Figure 2. Energy difference between the diquark 
condensate metastable state and the equilibrium 
state. 

order 1/g cannot remove such energy difference 
if g is sufficiently large. Hence, diquark conden- 
sation cannot take place in the strong coupling 
region. If, on the other hand, diquarks condense 
at some finite chemical potential in the weak cou- 
pling regime, there must be a critical value of 
the coupling, g^, such that diquark condensation 
takes place in some interval of the baryon den- 
sity for g < gc, but diquarks do not condense 
at any density for g > gc- This means that the 
interval of baryon densities (in lattice units) at 
which diquark condensation occurs as a stable 
state shrinks as g increases, and vanishes at gc- 
For g > gc the state with diquark condensate sur- 
vives as a metastable state. 
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